Fluctuation-Induced Interaction between Randomly Charged Dielectrics 
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Monopolar charge disorder effects are studied in the context of fluctuation-induced interactions 
between neutral dielectric slabs. It is shown that quenched bulk charge disorder gives rise to an 
additive contribution to the net interaction force which decays as the inverse distance between the 
slabs and may thus completely mask the standard Casimir-van der Waals force at large separations. 
By contrast, annealed (bulk or surface) charge disorder leads to a net interaction force whose large- 
distance behavior coincides with the universal Casimir force between perfect conductors, which 
scales as inverse cubic distance, and the dielectric properties enter only in subleading corrections. 
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Recent ultrahigh sensitivity experiments on Casimir 
(zero temperature and ideally polarizable surfaces) and 
van der Waals (finite temperature and non-ideally po- 
larizable surfaces) interactions between surfaces in vacuo 
[TJ 12] have highlighted the need for an accurate assess- 
ment of the possible electrostatic contribution to the to- 
tal interaction when the surfaces bear a disordered charge 
distribution [JT . The surface charge distribution can have 
various origins. In the so-called patch effect, the vari- 
ation of the local crystallographic axes of the exposed 
surface of a clean polycrystalline sample can lead to a 
variation of the local surface potential j4j- These varia- 
tions are of course sample specific and depend heavily on 
the method of preparation of the samples. The electro- 
static forces due to this surface potential disorder cannot 
be eliminated by grounding the two interacting surfaces. 
A similar type of surface charge disorder can also be ex- 
pected for amorphous films deposited on crystalline sub- 
strates. Surface annealing of these films can produce a 
grain structure of an extent that can be larger than the 
thickness of the deposited surface film [5 . In addition, 
adsorption of various contaminants can also influence the 
nature and type of the surface charge disorder. 

Here we assess the effect of various types of monopolar 
charge disorder on the interaction between two macro- 
scopic surfaces, delimiting two semi-infinite net-neutral 
dielectric slabs, separated by a layer of vacuum or an 
arbitrary dielectric material (Fig. [ij. Since the nature 
and distribution of the charge disorder in any of the ex- 
periments is in general not known, we consider different 
a priori models for the distribution of disorder. Specifi- 
cally, we assume that the charge disorder originates from 
randomly distributed monopolar charges which may be 



present both in the bulk and on the bounding surfaces 
and can be either annealed or quenched. It turns out 
that the type and the nature of the disorder has im- 
portant consequences for the total interaction between 
apposed bodies and can even dominate or give a con- 
tribution comparable to the underlying Casimir-van der 
Waals (vdW) effect. Our main goal is thus to investi- 
gate the interaction fingerprint of the charge disorder 
and to compare its contribution to the total interac- 
tion with the zero-frequency Casimir-vdW interaction 
between macroscopic surfaces. This may in turn help 
in assessing whether the experimentally observed inter- 
actions can be interpreted in terms of disorder-induced 
effects or pure Casimir-vdW interactions. 

We consider two semi-infinite slabs of dielectric con- 




FIG. 1: (Color online) A typical experimental set up (bottom) 
is modeled with a plane-parallel system (top) of two dielectric 
slabs (half-spaces) of dielectric constant Sp interacting across 
a medium of dielectric constant Sm- The charge distribution 
in the bulk of the slabs and on the two bounding surfaces at 
z = ±_D/2 has a disordered component (shown schematically 
by small light and dark patches) with zero mean but finite 
variance, and may be either quenched or annealed in nature. 
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stant Ep and temperature T with parallel planar inner 
surfaces (of infinite area S) located normal to the z axis 
at z = ztD/2, where D is thus the distance between their 
surfaces (see Fig. [l]). The inner gap is filled with a ma- 
terial of dielectric constant £„. We shall assume that 
the dielectric slabs have a disordered monopolar charge 
distribution, p(r), which may arise from charges residing 
on bounding surfaces [ps(r)] and/or in the bulk [p6(r)], 
i.e. p{r) = Ps{r) + Pb{'^)- The charge disorder will be 
taken to be of zero mean (i.e., the slabs are net neutral) 
and Gaussian-distributed ^ with no correlation in space 
i.e., ((p(r)p(r'))) = g{r)5{r — T') [where ((• • • )) denotes the 
disorder average] . The total correlation is the sum of the 
surface and bulk correlations g{v) = gsir) + gf,(r). For 
the slab geometry considered here, the charge distribu- 
tion is assumed to be statistically invariant in the plane 
of the dielectrics but with a variance dependent on z as 



gs{r) = gsel[S{z- 
5'b(r) = 



9bel 




-D/2) + S{z 

\z\ > D/2, 
\z\ < D/2, 



(1) 
(2) 



where eg is the elementary charge. It is worth mention- 
ing how this sort of disorder distribution might arise. If 
the bulk material has charge impurities at the sites 
distributed uniformly and independently with density Ub 
and charges Qi = ±q{,eo with equal probability, then we 
clearly have pb{r) = J^i 9i'5(r - ^i) and find ((pb(r))) = 
and {{pb{r) pb{r))) = qleQnb5{r — r'). We can thus make 
the correspondence gt = Qb^b- Similarly, one can make 
the correspondence gs = Qsns with Ug being the surface 
density of impurity charges ig^ on bounding surfaces. 

The partition function for the classical Casimir~vdW 
interaction (the zero-frequency Matsubara modes of the 
electromagnetic field) may be written as a functional in- 
tegral over the scalar field <j>{r), 



Z[pir)] ^ J [V<p{r)] 



-/35[0(r);p(r)] 



(3) 



with /3 — l/k^T and the effective action 

5[0(r);p(r)]= J dr [^eoe{r) {V^{r)f+i p{r)^{r)], (4) 

where e(r) = Ep for \z\ > D/2 and Em otherwise. In order 
to evaluate the averaged quantities such as the effective 
interaction, one needs to average the partition function 
over different realizations of the disordered charge dis- 
tribution, p{r) [7]. It is thus important to distinguish 
between quenched and annealed disorder that involve dif- 
ferent averaging schemes. For quenched disorder, the 
disorder charges are frozen and cannot fluctuate; for an- 
nealed disorder the charges can fluctuate and, in partic- 
ular, the charge distribution in the two slabs can adapt 
itself to minimize the free energy of the system. These 
two different disorder types lead to very different physical 
behaviors as will be demonstrated below. 



Let us first consider the quenched disorder model, 
where one must take the disorder average over the sample 
free energy, lnZ[/3(r)], in order to calculate the averaged 
quantities. The free energy of the quenched system. 



/3^qucnchcd = -((lnZ[p(r) 
may be calculated from Eq. (|3|) as 



(5) 



/?^qucnched= ^TrlnG-i(r,r')+^Tr{.9(r)G(r,r')}, (6) 



where G(r,r') is the Green's function defined via 
eoV- [£(r)VG(r,r')] - -(5(r-r'). 



(7) 



^ , we recognize 
vdW interaction 
The second term. 



In the first term of Eq. 
the standard zero-frequency 
/3.Fvdw = iTrlnG-i(r,r'). 
f3Tg = ^Tr{(7(r)G(r, r')}, represents the contribu- 
tion of the quenched charge disorder, which turns out to 
be additive in the free energy. 

The quenched expression ([6| is valid for any arbitrary 
disorder variance g{r). We now particularize to the case 
of planar dielectrics by using Eqs. ([T|) and in which 
case the zero-frequency vdW contribution per unit area. 



/3J-"vdw _ 1 
~ 2 
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In (1 - A 



2 ^-2QD 
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yields the standard vdW force, /, 



vdW 



/3/vdw 

s 



Li3(A^) 
87rL>3 ■ 



'dT^dw/dD, as 



(9) 



The dielectric jump parameter is defined as A = (Ep — 
£m) /(sp+^m) and Li3(-) is the trilogarithm function. The 
bulk and surface disorder contributions are obtained as 



2gblB£ni^ 
25s^Bem|ln(l 
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at all separations D with Ib — /3e§/(47r£o) — 56.8 nm 
being the Bjerrum length in vacuum at room temperature 
(T — 300 K). The quenched contribution from the bulk 
disorder (first term) in Eq. (10) is in principle infra-red 
divergent, however the corresponding force is finite. The 

total force, /quenched = -9J^quenched/9D, thuS foUowS aS 



/-^/quenched 

S 



gbln^ ( '^.gshe 



2EpD 



|ln(l - A2) 
+ ep)2A 



1 Li3(A2) 
(11) 

Here we obtain a sequence of scaling behaviors of differ- 
ent origins: a leading 1/D term due to the quenched bulk 
disorder, a subleading term from the surface charge 
disorder, and the pure vdW term that goes as and 
prevails in the absence of charge disorder. The disorder 
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FIG. 2: a) Ratio of the total force (111 to the zero-frequency vdW force ([ojl between net-neutral dielectric half-spaces (in 
vacuum, e™, = 1) bearing quenched monopolar charge disorder for fixed bulk and surface disorder variances gt — 5x 10~* nm~^, 
ffa — ^iid different dielectric constants £p = 2, 5, 10, 20, 40, 80, 100 (dashed curves from top), b) Magnitude of the rescaled 
total force (111 in the quenched case as a function of the rescaled distance for fixed Sp = 10, Qs = and various bulk disorder 
variances (/(, = 10~^, lO"'^, 10~*, 10~^, 10~^°, lO"'^^ nm~^ (dashed curves from top). Solid curve is the pure vdW force Inset 
is the ratio of the total force to the vdW force ^ for the same range of D. c) is the same as (b) but for annealed disorder 
(top panel, from Eq. ( |16[ l). Annealed curves stay close to one another and are bracketed by the perfect conductor result (Eq. 
(20l, top solid line, labeled by C(3)) for large disorder and the vdW result for no disorder (Eq. bottom solid line, labeled 
by Li3(A^)) as seen more clearly from the force ratio sho-wn in the bottom panel, (b) and (c) are plotted in log-log scale. 



contributions (first and second terms in Eq. (Ill) are 



attractive "when the dielectric mismatch A > (e.g., for 
t-wo dielectrics slabs in vacuum) and repulsive other"wise 
(e.g. , for the t"wo surfaces of a single slab in vacuum) . One 
might expect that globally electroneutral slabs would ex- 
hibit a dipolar-like interaction force on the leading or- 
der rather than the monopolar forms 1/D (or l/D^) ob- 
tained for the bulk (or surface) charge distribution. The 
physics involved is indeed subtle as the disorder terms 
result from the self-interaction of the charges -with their 
images (-which foUo-ws from G(r,r), Eq. ([6]), and only 
in a dielectrically inhomogeneous system) and not from 
dipolar interactions ("which come from an expansion of 
G{r,r') "when |r — r'| is large). Statistically speaking each 
charge on average (as any other charge has an equal prob- 
ability of being of the same or opposite sign) only sees 
its image, thus explaining the leading monopolar form in 
the net force. 

The remarkable result is however the relative impor- 
tance of the disorder-induced forces, which exhibit a 
much weaker decay with the separation, D. They may 
thus completely mask the standard Casimir~vdW force at 
sufficiently large separations depending on the dielectric 
constants and the disorder variances as shown in Figs.|2|i 
and ^p; the bulk disorder variance is chosen here within 
the typical range gt, ~ 10"-'^^ — 10~^ nm~^ (corresponding 
to impurity charge densities of 10^° — 10^^ eo/cm'^) [8 . 

For relatively small surface disorder {gs <C 56 ?b), the 
anomalous l/D behavior is predicted to dominate the 
vdW behavior beyond the crossover distance 



£pLi3(A^) 
'ingtlBA 



1/2 



(12) 



which, for typical parameter values, covers the experi- 
mentally relevant range of distances from a few hundreds 
of nm to several microns (Fig. [2Jd, main set). For strong 
surface disorder (g^ S> 56?b), on the other hand, one ex- 
pects the 1 /D^ behavior to dominate beyond 

(£„ + £p)'^ALi3(A") 

' 16^5s;Be™|ln(l-A2)|- ^ ^ 

So far we have only examined the effects from the 
quenched disorder. In reality one may encounter disor- 
dered charges with some degree of annealing. A possible 
idealization is to assume that the disorder is completely 
annealed (the intermediate partially annealed cases are 
also analytically tractable as shown recently ^ but will 
not be considered here). In the annealed disorder model, 
the disorder average is taken over the sample partition 
function, Z[/9(r)]; hence, the free energy of the system. 



/3.Fanncalcd = - In ((Z[p(r)])) , 



(14) 



may be evaluated as 
/3^an„eaicd = ^Trlu [G-\r,r') + l3gir)Sir ~ r')] . (15) 

Note that, unlike the quenched result in Eq. ([6]), the 
disorder and the pure Casimir-vdW contributions can 
not be separated in general when the disorder is annealed. 

In the case of two interacting slabs with the surface 
and bulk disorder variances ([TJ and (|2|, the modified 
inverse Green's function G'~^(r, r')+f3g{r)5{r—r') may be 
evaluated explicitly and the fluctuational trace-log term 
may be calculated by the standard methods [TU] as 



annealed 



/|^ln(l-A^(g)e-n (16) 
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with 

Ag(g) 



SmQ -ep^/Q^ + 4:TTlBgb/£p - 47r^Bg^ 
_^ inlsgb/sp + inlsgs 



(17) 



Let us first consider the large-distance behavior of the 
net annealed force. For strong annealed bulk disorder 
(<?s <C ffb^B), we obtain the asymptotic behavior 



annealed 



C(3) , / 3C(3)e„ 



1 



(18) 



which is expected to hold for D » iSm/ ^T^gb^B^p [i-e., 
-D ^ 70 nm for gb = 10^® nm~'^ and Ep = 10 in vacuum]. 
While for weak bulk disorder (g^ ^ Sb^e), we obtain 



annealed 



C(3) / 3C(3) 



1 



87ri:>3 yiQ-K^gslBj D^' 



(19) 



which is expected to hold for D ^ Sem/i^ngslB) [i-e., 
-D 3> 80 nm for gs — 10^* nm^^ in vacuum]. Obvi- 
ously, material properties disappear in the leading-order 
total force between arbitrary dielectrics bearing annealed 
charge disorder and one ends up with the universal at- 
traction as one would expect for two perfect conductors 



PL 



annealed 



C(3) 
87rL>3 



D 



(20) 



These asymptotic behaviors are also obtained for strong 
disorder (gb or gs oo). On the contrary, for weak dis- 
order (gb and gs — *■ 0) or for vanishing separation, one re- 
covers the non-universal vdW force ([9| as the asymptotic 
behavior. It is thus interesting to note that the force in 
the annealed case is bounded between these two limiting 
results, i.e. Eqs. ^ and (|20| (Fig. [2j;, solid lines). The 
above results demonstrate the intuitive fact that dielec- 
tric slabs with annealed charges tend to behave asymp- 
totically as perfect conductors and, unlike the quenched 
case, the algebraic decay of the leading contribution to 
the net force remains unchanged. The deviations due to 
material properties and the disorder variance contribute 



a repulsive subleading force (Eqs. (18) and (19l). 

For the experimental sphere-plane geometry [2\ a naive 
application of the proximity force approximation [1 to 
the results derived above would lead to forces with the 
leading behavior '--^ InD -I- 0{D~^) in the quenched case 
(from Eq. (11 1) and ^ + 0{D~^) in the annealed 
case (from Eqs. (18 1 and (19l). Thus, an effective 



scaling exponent (defined as D"") of a < 1 (consis- 
tent with recent experimental observation of a resid- 
ual electrostatic force scaling as D~'^-^ [TT]) may be ob- 
tained in the quenched case, both with the bulk disorder 
(plane-plane and sphere-plane geometry) and the sur- 
face disorder model (sphere-plane geometry). A more 
detailed comparison with force measurements should be 
attempted once the experimental and methodological un- 
certainties surrounding experiments are sorted out (see 
Ref. pi and published comments). 



In conclusion, we have studied the influence of charge 
disorder on the fluctuation-induced interaction between 
net-neutral dielectric slabs bearing random quenched 
and/or annealed monopolar charges on their bounding 
surfaces and/or in the bulk and compared it with the 
zero- frequency Casimir-vdW interaction between them. 
Quenched disorder leads to an additive contribution to 
the total interaction force that scales as 1/D (or 1/-D^) 
for bulk (or surface) charge disorder, may be attractive or 
repulsive and depends on the dielectric constants of the 
materials. By contrast, annealed disorder leads to a net 
attractive force, which is universal and decays as 
at large separations. Thus, the main flngerprint of the 
annealed disorder is that the total force remains intact 
in this case as the dielectric constants are varied. This 
could help distinguish this type of interaction from the 
pure Casimir-vdW interaction between dielectrics with 
no disorder as well as from the interaction due to the 
quenched disorder, Eq. (111. These two latter cases can 



in turn be distinguished by monitoring the dependence 
on the separation of the net interaction, which for the 
quenched disorder exhibits a much weaker decay. Note 
that the disorder effects are compared here with the zero- 
frequency Casimir-vdW interaction. The precise correc- 
tion presented by the higher-order Matsubara frequencies 
is very material specific, but its magnitude (relative to 
the zero-frequency term) is typically small for the most 
part of the separation range considered here and remains 
negligible in comparison with the quenched terms in (11). 

We thank M. Kardar and W.J. Kim for useful dis- 
cussions. This research was supported in part by the 
NSF (Grant No. PHY05-51164). D.S.D. acknowledges 
support from the Institut Universitaire de France. R.P. 
acknowledges support from ARRS. A.N. is a Newton In- 
ternational Fellow. 



[1 
[2 

[3 
[4 
[5 

[6: 

[7 

[s: 

[9 

[lo: 



M. Bordag, U. Mohideen and V.M. Mostepanenko, Phys. 
Rep. 353, 1 (2001). 

W.J. Kim et al, Phys. Rev. A 78, 020101(R) (2008); 
79, 026102, (2009); Phys. Rev. Lett. 103, 060401 (2009); 
R.S. Decca et al, Phys. Rev. A 79, 026101 (2009); S. de 
Man, K. Heeck, and D. lannuzzi, ihid 79, 024102 (2009). 
C.C. Speake and C. Trenkel, Phys. Rev. Lett. 90, 160403 
(2003). 

L.F. Zagonel et al., Surface and Interface Analysis 40, 
1709 (2008). 

Z.H. Liu, N.M.D. Brown and A. McKinley, J. Phys.: 
Condens. Matter 9, 59 (1997). 

The Gaussianity assumption is unimportant for quenched 
charges and only the correlation function ((p(r)p(r'))) 
needs to be given. 

A. Naji, R. Podgornik, Phys. Rev. E 72, 041402 (2005); 
R. Podgornik, A. Naji, Europhys. Lett. 74, 712 (2006). 
K.C. Kao, Dielectric Phenomena in Solids (Elsevier Aca- 
demic Press, San Diego, 2004); L.P. Pitaevskii, Phys. 
Rev. Lett. 101, 163202 (2008). 

Y.S. Mamasakhlisov, A. Naji and R. Podgornik, J. Stat. 
Phys. 133, 659 (2008). 

R. Podgornik, J. Chem. Phys. 91, 5840 (1989). 



[11] W.J. Kim et aZ., |arXiv:0905.3421|^l (2009) 



